Radiative Lifetime of Excitons in Carbon Nanotubes 
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We calculate the radiative lifetime and energy bandstructure of excitons in semiconducting carbon 
nanotubes, within a tight-binding approach. In the limit of rapid interband thermalization, the 
radiative decay rate is maximized at intermediate temperatures, decreasing at low temperature 
because the lowest-energy excitons are optically forbidden. The intrinsic phonons cannot scatter 
excitons between optically active and forbidden bands, so sample-dependent extrinsic effects that 
break the symmetries can play a central role. We calculate the diameter-dependent energy splittings 
between singlet and triplet excitons of different symmetries, and the resulting dependence of radiative 
lifetime on temperature and tube diameter. 
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The potential usefulness of carbon nanotubes as an op- 
tical material depends sensitively on the luminescence ef- 
ficiency, which in turn depends on the radiative lifetime. 
Yet while the opticalproperties of nanotubes have been 
studied intensely Q, H H H, |Ili,|7lli|i|il|ll[ll[l3, 
little is known about the radiative lifetime, and its depen- 
dence on temperature and tube diameter. It is by now 
well recognized that optical absorption and emission in 
nanotubes are dominated by excitons 0, 0, 0| , posing 
an additional challenge. 

Here we calculate the exciton bandstructure and ra- 
diative lifetime in carbon nanotubes. When interband 
thermalization is efficient, we find that the lifetime has 
a minimum at intermediate temperatures, with optical 
recombination becoming slower both at very low and 
very high temperature. This unusual behavior occurs be- 
cause the lowest-energy excitons are optically forbidden 
by symmetry 0, Q| . We report the energy differences 
between singlet and triplet excitons of different symme- 
tries, for a range of nanotube diameters. We also present 
the resulting temperature dependences of the radiative 
lifetime. Both the energies and temperature dependences 
are well described by simple scaling formulas. The opti- 
cally allowed band has a non-analytic dispersion due to 
exchange, while the other bands have a mass enhance- 
ment of ~50% over free electron-hole pairs. 

We find that symmetry plays a special role in this sys- 
tem. For photoluminescence of an ideal nanotube (and 
neglecting spin-orbit coupling), phonons cannot scatter 
excitons into the optically forbidden low-energy states, 
so optical emission is very efficient. However, extrinsic 
effects that break the ideal symmetry can qualitatively 
change the behavior, especially at low temperature. Such 
effects can also lead to multiple peaks in emission spec- 
tra with characteristic temperature dependence, perhaps 
explaining experimental observations 0, 0| . 

The quantum efficiency 77 of light emission corresponds 
to the fraction of excitons decaying by radiative channels: 
V — Tnr/{Tr + T„r), whcrc Tr and T„r are the radiative 
and nonradiativc lifetimes. The total decay rate (i.e. in- 



verse lifetinie) is_n~( — +'''nr ■ It been measured 
[itI fl^ Hfll EH Elf , and environmental effects have been 
recently published But only one measurement of 

radiative lifetime has been reported [20j . 

A photoexcited state has total momentum q = 
g(photon), hereafter approximated as g = 0. It is rapidly 
scattered by phonons into a different state with g 7^ 
01, and cannot decay directly. Therefore what is rele- 
vant is not the lifetime of a specific excitonic state, but 
rather the lifetime of an exciton as it scatters from one 
state to another. If thermalization is inefficient, then the 
lifetime depends on the details of the excitation process 
and subsequent scattering pathways. 

We therefore focus first on the case where the exciton 
thermalizes on a timescale much shorter than the decay 
time [13 III III in E3- Then at low excitation den- 
sity there are well-defined (thermally averaged) radiative 
and non-radiative lifetimes. The radiative decay rate (i.e. 
inverse radiative lifetime) is 23j 
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and a similar formula applies for nonradiativc decay. 
The index i' = {q,j,S,a) may include summation over 
exciton wavevector q, angular momentum j (i.e. dis- 
crete wavevector around the tube), spin index a (la- 
belling singlet and triplet states), and exciton state S 
(labelling the bound and continuum states for a given 
q and a); or it may be restricted to those states ac- 
cessible by symmetry-allowed scattering mechanisms, as 
discussed below. The partition function is Z(T) = 
J2iyGxp{—E^/kBT), summed over the same range of v. 
The radiative decay rate for each excitonic state is 
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where n is the refractive index (hereafter n — 1.4). We 
calculate the exciton energies Eu and the correspond- 
ing oscillator strengths fjy using a tight-binding Hamilto- 
nian E^ with t — 3.0 eV, and solving the Bethe-Salpeter 
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FIG. 1: Four lowest-energy exciton bands of a (19,0) nan- 
otube {d — 1.5 nm) embedded in medium having e = 2 [13 . 
for (a) singlet and (c) triplet spin states. The black solid 
curves are doubly degenerate exciton bands having nonzero 
angular momentum. The blue dashed curves are even par- 
ity (dipole forbidden) exciton bands having zero angular mo- 
mentum. The red dotted curves are the odd parity (dipole 
allowed) exciton bands having zero angular momentum. The 
density of states for (b) singlet and (d) triplet states are shown 
by the black curves. The red dotted curve in (b) shows the 
partial density of states for dipole allowed singlet exciton. The 
change in scale of 1/3 in (d) is for triplet degeneracy. 



Equation (BSE) 113, as in Refs. [T1[T3|, taking the sur- 
rounding medium to have dielectric constant e = 2. 

In an ideahzed picture, one exciton can capture much 
of the spectral weight of the entire nanotubes . But in 
reality, phonon scattering limits the exciton coherence to 
a finite length scale, and spreads the oscillator strength 
of the exciton over a finite linewidth |23|. As long as 
the linewidth is substantially less than fc^T 0|, we ob- 
tain good accuracy from Eq. without including this 
linewidth explicitly. 

For a single ID exciton band with parabolic disper- 
sion, the temperature dependence of the radiative life- 
time [Eq. (^] is Tr oc T^l"^ Q^. However, as shown in 
Fig. n here there are several low- lying bands , with the 
lowest being optically forbidden |l2l [ij] . Thus it is im- 
portant to understand the exciton energies and symme- 
tries in some detail. 

The conduction and valence bands in carbon nan- 
otubes are doubly degenerate. For free electron-hole 
pairs, this gives four distinct but degenerate pair exci- 
tations. In excitons, the Coulomb interaction partially 
lifts the fourfold degeneracy. The dispersions of the four 
lowest-energy singlet excitons are shown in Fig. The 
lowest-energy exciton at 5 = has zero angular mo- 
mentum and even parity, so it is dipole forbidden [l2l |. 



The second exciton also has zero angular momentum, 
but with odd parity, and it gives rise to the E\\ optical 
transition. The next two excitons at g = lie above the 
E\\ transition energy. They are degenerate, with posi- 
tive and negative angular momentum, and are therefore 
optically forbidden; but they can be observed in phonon- 
assisted absorption [l^ . 

The optically forbidden excitons all have similar effec- 
tive masses, m* ^ 1.5(me + m;i), where me = m;i is 
the effective mass of an electron or hole. The 50% mass 
enhancement is roughly independent of tube diameter. 

In contrast, the dispersion of the optically active ex- 
citon is extremely steep and non-parabolic, going as 
(?^ln|g|. A similar logarithmic contribution to the elec- 
tron self energy was found in 2D graphene [2^. This 
anomalous dispersion occurs because of the exchange en- 
ergy, which for the other bands is either zero by symme- 
try, or is a smooth function of q. 

In photoluminescence experiments, excitons are opti- 
cally pumped into a singlet spin state of odd parity. For 
an ideal nanotube (and neglecting spin-orbit coupling), 
we find that phonons cannot scatter an exciton between 
states of different parity or different spin. Therefore the 
exciton will thermalize within the odd-parity spin-singlet 
bands (sohd black and dotted red in Fig. ^), and the 
sums in Eq. and Z(T) need to be restricted accord- 
ingly. The resulting behavior is shown in Fig. At 
low temperature, the behavior is well approximated as 
Tr cx shghtly different than T^/^ dependence [lil 

due to the nonparabolic dispersion. At higher temper- 
ature, the behavior is quite different, reflecting the con- 
tribution of the higher-energy bands to Z{T); it is well 
approximated by oc T"^/^. 

This ideal picture, while appealing in its simplicity, 
may not be directly applicable to experiment. There are 
many effects that tend to break the symmetries of the 
perfect nanotube: coupling to the environment, struc- 
tural distortions, finite size and associated end effects, 
and the presence of defects or impurities. These effects 
allow scattering and thermalization between bands of dif- 
ferent parity and/or different spin. This tends to sup- 
press luminescence at low temperature, because excitons 
thermalize to a low-energy band that is optically forbid- 
den. At the same time, luminescence is suppressed at 
high temperature because a larger range of q becomes 
thermally accessible. These effects enter via Z{T) in 
Eq. id. 

Fig. [2b shows the case of scattering and thermaliza- 
tion between different parity, but not between different 
spin states; and Fig. ^ shows the case of thermaliza- 
tion among states of all parities and spins. We find that 
the temperature dependence of Tr is very similar in all 
cases where thermalization includes a low-energy forbid- 
den band, except for an overall scaling of t and T. This 
is shown in Fig. [2^, where all the curves in Figs, ^p-c 
collapse nearly to a single scaled curve. A simple in- 
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FIG. 2; Exciton radiative decay rate t^^ [Eq. as a func- 
tion of temperature, in a medium having e = 2 Curves in 
(a-c) are for zig-zag tubes with diameters 1.0 nm (dotted red), 
1.5 nm (soHd blue), and 2.0 nm (dashed green). Fine dotted 
black hne in (a) shows fit of 1.0 nm curve to T-^-^"" at low T 
and at high T. In (a) only singlet states of odd par- 

ity are included; (b) assumes thermalization between singlet 
states of even and odd parity; and (c) assumes thermahzation 
between both singlet and triplet states of both parities. Ef- 
fect of symmetry mixing is illustrated in (d) for 1.5 nm tube 
as in (b), but including a 5% mixing of spectral weight into 
the even-parity band. Solid blue curve shows decay rate from 
main _Eii peak, which is only very slightly reduced relative 
to (b). Dotted green line is contribution from nominally for- 
bidden exciton 3.5 meV lower in energy. Dashed red line is 
total radiative decay rate. Temperature dependence of dotted 
green curve depends on g-dependence assumed for mixing, (e) 
Scaled radiative lifetimes for all six curves from (b) and (c). 
The empirical scaling formula Eq. Q is shown by a dotted 
black-on-white curve. 

terpolation formula can describe this scaled temperature 
dependence: 

— - TfT- CXp . (3) 

' min min \ / 

The diameter dependence of Tmin and Tmin are shown 
in Fig.|SJi and|2b respectively, along with the best fits to 
the simple empirical relations: 

T"min « To + ad, Tmin ~ ^- (4) 

For small-diamctcr (large-bandgap) tubes, the maximum 
radiative decay rate t^I^ is higher, and this maximum 
rate occurs closer to room temperature. 
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FIG. 3: Diameter dependence in zig-zag semiconducting nan- 
otubes of (a) the minimum radiative lifetime, (b) the corre- 
sponding temperature, and (c) the energy splitting between 
the optically active exciton and the lowest energy exciton, 
along with best fits to Eq. Black circles are the results, 

when both triplet and singlet states are included; red squares 
are with only singlet excitons included. Thus the black cir- 
cles in (c) correspond to the singlet-triplet splitting of the 
dipole-allowed states at g = 0. Solid curves are interpolation 
formulas, Eq. Q and Eq. @. Note the xlO multiplication 
factor for the red squares in (a). 

Tmin is approximately proportional to the energy differ- 
ence between the optically active exciton and the 
lowest-energy exciton that participates in thermalization 
(Edo)- This difference En — E^o is shown in Fig.|3;. The 
singlet-triplet exchange splitting of the dipole-allowed ex- 
citon, as well as the even-odd parity splitting of the sin- 
glet exciton, are both inversely proportional to the tube 
diameter: 

En-Edo-^. (5) 

The singlet-triplet splitting is expected to scale with the 
size of the exciton, which is proportional to the tube di- 
ameter The splitting between the even and odd 
parity excitons is a factor of five smaller than the singlet- 
triplet exchange splitting. 

In Fig. 131 Tniin, Tniin, and Ell — Edo all exhibit alter- 
nations in the diameter dependence, due to the system- 
atic difference between tubes having mod(n-m,3)=l or 
2. While only zigzag tubes are shown here, we find very 
similar diameter dependence for other chiralities. 

Equations provide an estimate of the tempera- 

ture dependent radiative lifetime in this tube diameter 
range, whenever thermalization with optically forbidden 
symmetries occurs. The best-fit values of the three pa- 
rameters are: tq = 12 ns, a = 24 ns/nm, and (3 = 272 
K-nm, if both singlet and triplet excitons contribute to 
Z{T); and tq = ns, a = 2.2 ns/nm, and (3 — 99 K-nm 
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if only singlet excitons contribute. The singlet-triplet 
exchange splitting and the even-odd parity splitting at 
q = are given by Eq. jSJl, with 7 = 27 meV-nm and 
7 = 5.5 meV-nm respectively. 

At sufficiently low temperature, additional complica- 
tions arise. Thermalization by acoustic phonons becomes 
less effective [s^, so it is less likely that this equilib- 
rium picture will apply. Also, even weak disorder from 
environmental interactions would lead to localization of 
band-edge excitons. Because localized excitons that are 
close in energy can be well separated spatially, they typ- 
ically are not in thermal equilibrium at low T. This will 
give a finite linewidth determined by the strength of the 
disorder, even at T = 0. 

Interband scattering necessarily implies some mixing of 
the bands; so in the cases of inter-symmetry scattering 
(Fig. |23-c), the optically forbidden bands are no longer 
strictly forbidden. Rather, there is some small but finite 
radiative decay rate associated with every participat- 
ing band. For thermalization at sufficiently low T, only 
the lowest-energy of these bands will have significant oc- 
cupancy. Then the emission from this "nearly-forbidden" 
band, however small, will dominate over thermal excita- 
tion into the optically allowed band. 

In this case, the emission spectrum (not shown) has 
two peaks, one from the optically allowed band, and 
one at lower energy from the nominally forbidden band. 
The latter may be negligibly weak at high temperature, 
but becomes dominant at sufhciently low temperature 
(as long as thermalization occurs before decay). We il- 
lustrate this in Fig. [21, for one particular model of the 
mixing. 

Thus the behavior is sensitive to deviations from ide- 
ality, which will in general vary among different experi- 
ments. Single-tube spectroscopy at low temperature pro- 
vides a powerful diagnostic of such extrinsic effects. Any 
distinct emission peak at lower energy than the main 
peak implies some symmetry breaking. If the peaks do 
not correspond to the few distinct bands of Fig. ^ they 
are likely to correspond to excitons that are localized 
(e.g. at defects or deformations). A finite peak width as 
T — » suggests disorder (along with the expected break- 
down of thermalization) . Recent measurements are sug- 
gestive in showing multiple peaks, with spectral weight 
shifting from hig h- to low-energy peaks with decreasing 
temperature jl5lllf]l|. 
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